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We obtain the in-medium effective potential of the three-flavor Polyakov-quark-meson model as 
a real function of real variables in the Polyakov loop variable, to allow for the study of all possible 
minima of the model. At finite quark chemical potential, the real and imaginary parts of the effective 
potential, in terms of the Polyakov loop variables, are made apparent, showing explicitly the fermion 
sign problem of the theory. The phase diagram and other equilibrium observables, obtained from 
the real part of the effective potential, are calculated in the mean-field approximation. The obtained 
results are compared to those found with the so-called saddle-point approach. Our procedure also 
allows the calculation of the surface tension between the chirally broken and confined phase, and the 
chirally restored and deconfined phase. The values of surface tension we find for low temperatures 
are very close to the ones recently found for two-flavor chiral models. Some consequences of our 
results for the early Universe, for heavy-ion collisions, and for proto-neutron stars are discussed. 

PACS numbers: 25.75.Nq, 64.60.Q-, 11.30.Rd, 98.80.Cq 



I. INTRODUCTION 

The determination of the properties of strongly inter- 
acting matter at high temperatures and baryon densities 
is one of the main goals of today's high-energy physics. 
However, not many tools are available for such a difficult 
task. Due to the fermion sign problem of quantum chro- 
modynamics (QCD) at nonvanishing (real) quark chemi- 
cal potential, Monte-Carlo calculations on the lattice are 
not feasible in this regime due to the lack of an impor- 
tance sampling procedure that is free of ambiguities [HQ. 
(Nevertheless, the authors of [|[ have been able to inves- 
tigate the curvature of the T — fi phase diagram at low 
values of chemical potential on the lattice, with physical 
quark masses and in the continuum extrapolation, indi- 
cating that the crossover temperature decreases with the 
chemical potential for /Lt/T c <C 1.) An alternative ap- 
proach is that of chiral effective models for QCD, which 
have been successfully utilized for many decades [1, [j| . 
In this work, we shall adopt one such model, the three- 
flavor Polyakov-quark-meson (PQM) model [6^. It has 
become quite popular in the last years due to its close 
relationship with the linear sigma model and its agree- 
ment with results from lattice calculations of thermo- 
dynamical quantities at zero baryon chemical potential 
[Tol [Tlj . most importantly the crossover nature of the 
quark-hadron transition at vanishing chemical potential 
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As discussed for example in [13l4l5| . not only QCD, 
but also effective models that possess gauge degrees of 
freedom do present the sign problem at finite chemical 
potential /x, even at the mean-field level. In such mod- 
els, this is manifest in the appearance of an imaginary 
part of the in-medium effective potential in equilibrium 
at fj, 0. Since the resulting effective potential at fi- 
nite chemical potential is a complex function of complex 
variables, special care must be taken with respect to the 
meaning of a minimization procedure that leads to the 
state of equilibrium of the system at any given nonzero 
temperature T and chemical potential /x. Following an 
approach similar to the one in [l5j (in the context of the 
PNJL model), we propose a change of variables, followed 
by a simple approximation in the PQM model that ren- 
ders the in-medium effective potential a real function of 
real variables. As a consequence, the effective potential 
is, in this approach, a real function that possesses min- 
ima, as demanded by general field-theoretical arguments 
applied to systems in equilibrium. 

The implementation of such an approximation scheme 
(that simply accounts for neglecting the imaginary part 
of the effective potential) of course leads to differences 
in the predictions made, for example, with the so-called 
saddle-point approach (see, e.g., [16[). In this approach, 
the state of thcrmodynamical equilibrium is found by 
first restricting the Polyakov loop variables to real vari- 
ables. Second, all the extrema of the effective potential 
(i.e., points in which the derivatives of the effective po- 
tential with respect to all order-parameters vanish) are 
determined. At a third step, the state of equilibrium is 
chosen among the extrema as the one with the lowest 
value of the effective potential. It often happens that the 
chosen point is not a minimum, but a saddle-point of the 
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effective potential. In fact, in this approach, the effective 
potential often has no minima, but only saddle-points. 

In this work, we show how the differences between the 
two approaches are manifest in the T — fx phase diagram 
of the model. We also briefly discuss differences and simi- 
larities between several parametrizations of the Polyakov 
loop potential found in the literature (see, e.g., 17-[l9j). 
The saddle-point approach is particularly troublesome 
for the computation of the surface tension in the region of 
the phase diagram of the model where a first-order phase 
transition exists. The computation of the surface tension, 
and thus of nucleation rates, requires a precise definition 
and location of the minima of the model. This adds an- 
other important aspect to our suggested approach, which 
is that it consistently allows the evaluation of the surface 
tension between two phases of the model, which are sep- 
arated by a first-order phase transition. As discussed 
for example in (2p| - [29| . the surface tension is a crucial 
parameter for the dynamics of a first-order phase transi- 
tion in many scenarios of high-energy physics. Roughly 
speaking, if the surface tension between the stable and 
the metastablc phases in a given phase transition is too 
high, the phase conversion may be exceedingly slow. In 
practical terms, a high surface tension can dynamically 
suppress a first-order phase transition that would be oth- 
erwise allowed if only the bulk equilibrium thermody- 
namics were considered. For this reason, the evaluation 
of the surface tension from effective models of QCD is 
a relevant question about strongly interacting systems 
that can undergo a first-order phase transition, such as 
the ones produced in relativistic heavy-ion collisions at 
RHIC, NIC A and FAIR, and also in the interior of com- 
pact stars. This evaluation, however, is not available 
from first-principles QCD but can be performed in a rel- 
atively straightforward manner using an effective model 
for QCD (with or without Polyakov loop degrees of free- 
dom), such as the two-flavor quark- meson model [30j . 
the Nambu-Jona Lasinio (NJL) model [3l[, or the PQM 
model, which we will be studying in this work. Another 
goal in this work is to produce a result that at least over- 
estimates the surface tension in the three-flavor PQM 
model along the line of first-order phase transition in the 
T— fx plane. This is an interesting result by itself, but it is 
also an interesting example of a calculation of the surface 
tension in a theory with multiple order-parameters. No- 
tice that the calculation of the surface tension explicitly 
demands not only the localization of the minima of the 
effective potential, but also the precise form of the poten- 
tial between them. We argue that having a real effective 
potential as a function of real order-parameters is crucial 
for the calculation of the surface tension, as well as for 
any other quantity in full thcrmodynamical equilibrium. 

The structure of the paper is as follows. After a brief 
summary of the PQM model, in Sec. [Til we will discuss 
the necessary conditions on the effective potential for the 
definiteness of thcrmodynamical equilibrium states and 
how these apply to the model. We also show how the 
sign problem explicitly appears in the PQM model at 



the mean-field level and propose an approximation in or- 
der to circumvent it. In Sec. IIII[ for completeness, we 
start the section by reviewing the general computation 
of surface tensions in a first-order phase transition. Af- 
ter that, we discuss how to obtain an overestimate of the 
surface tension of a phase interface in the PQM model. In 
Sec. II VI wc show our results for the phase diagram for the 
PQM model and the surface tension. The implications 
of these results arc also discussed. Finally, in Sec. [V] wc 
give our conclusions and final discussions of the obtained 
results. 



II. THE POLYAKOV-QUARK-MESON MODEL 
AT FINITE T AND /i 

In order to incorporate aspects from the physics of 
chiral symmetry breaking and restoration, as well as 
from the confinement-deconfincmcnt phase transitions, 
we adopt an effective model that captures these main fea- 
tures of quantum chromodynamics. In the last years, two 
such models have gained much popularity, the Polyakov- 
loop extended Nambu-Jona Lasinio (PNJL) model fl6l . 
[l8l.l32[ [33j and the Polyakov-loop-extended quark-meson 
(PQM) model H^. These models may be considered as 
extensions of the well-known Nambu-Jona Lasinio model 
[334111 and the linear sigma model [H which 
provide an effective realization of the chiral symmetry 
breaking pattern. The physics of confinement is expected 
to be taken into account by coupling the quark and the 
meson degrees of freedom with the expectation value of 
the Polyakov loop. 

In this work, our choice will be a PQM model with 
2 + 1 quark flavors, as in Refs. @, 



A. The Polyakov loop 

Before discussing the model, we briefly introduce the 
order parameter for (de)confincment with which the sign 
problem enters in Polyakov-loop extended models. 

The Polyakov loop operator is a Wilson loop in tem- 
poral direction, 



7 = P exp 




dx A (a? ) 



(2.1) 



where P denotes path ordering, f3 = 1/T is the inverse 
of the temperature and Aq is the temporal component of 
the gauge field (42j. With a gauge that ensures the 
time independence of Aq, we can perform the integration 
trivially and the path ordering becomes irrelevant [43l . 
[HI, so that y = exp (z/L4o). In this form, it is trivial to 
see that the Polyakov loop variable, 
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is a complex scalar field $ = a + i/3. Furthermore, we 
can rotate the gauge field in the cartan subalgebra = 

^o 3 ^3+^o 8 ^8 [HI- Within this diagonal representation, 
we see that the adjoint Polyakov loop variable becomes 
simply $ = a — 

Under gauge transformations concerning center sym- 
metry, the Polyakov loop operator and its variable are 
multiplied with a center element Z, $ — > In pure 

SU(3) gauge theory, that can be considered as QCD 
with infinitely heavy quark masses, the confining phase 
is center-symmetric and, therefore, ($) = 0, while decon- 
finement is characterized by a finite value of the Polyakov 
loop expectation value, since center symmetry gets bro- 
ken spontaneously Real QCD with physical quark 
masses adds an explicit symmetry breaking term. 



B. Formulation of the Model 

Let us briefly recall a few aspects of the PQM model. 
The interested reader may refer to e.g. Refs. ffl. I9L l38l. l4l| 
for more details. 

The starting point for our analysis of the equilibrium 
states of the system is the in-medium effective poten- 
tial of the theorjQ as a function of the relevant order- 
parameters and thermodynamical control parameters. 
This effective potential is derived within a mean-field 
approximation to a theory with constituent quarks min- 
imally coupled to gauge fields and coupled to mesons 
via a Yukawa-type term. In the leading in-medium con- 
tribution from the fcrmions, the coupling to the color 
fields effectively becomes a coupling to the Polyakov loop 
field. The coupling with the mesons is translated into 
the masses of the quarks via spontaneous and explicit 
chiral symmetry breaking terms in the self-interaction 
Lagrangian density for the mesons. The self-interaction 
in the gauge sector is modeled by a potential energy for 
the Polyakov loop variable. 

In the chiral sector of the theory, the natural choice 
of order parameters for the Nf = 2 + 1 pattern of chiral 
symmetry breaking are the non-strange (cr x ) and strange 
(a y ) chiral condensates [4l|. In the Polyakov- loop sec- 
tor, the natural variables are the Polyakov loop itself, 
and its conjugate, 3>. With these choices, the effective 
potential is given by the sum of three terms, 

a = u(a x , g v ) + u($, i) + n q9 (a x ,a y ,&, $) (2.3) 

where the dependence on the thermodynamical control 
parameters (in this case the temperature T and the quark 



1 As long as wc are interested in states of thermal equilibrium of 
the system, which are described by homogeneous field configura- 
tions, the in-medium effective action is reduced to an in- medium 
effective potential. In the next section, where we study non- 
homogeneous states, we will have to consider the effective action 
instead. 



TABLE I. Values of constants to which the parameters of 
the mesonic potential are adjusted, according to Ref. and 
value of the constituent quark mass of the light (up and down) 
quarks that we use to fix the quark meson Yukawa coupling 
in Eq. l|2T0]) . 



Constant 


A 


Jk 


m,r uik rn v m v i 


m a 


mi 


Value [MeV] 


92 


110 


138 495 548 958 


400 - 550 


300 



chemical potential fx) is implicit in the last two terms of 
Eq. (|2.3p . Below, we describe in detail each of the con- 
tributions appearing in the thermodynamical potential 
Eq. (23). 

The first term of the thermodynamical potential (|2.3[) 
is the tree-level contribution from the mesonic degrees of 
freedom and reads 



U (a Xl <7 y ) 



+1 (2Ai + A 2 ) o x + I (2Ax + 2A 2 X 



2V2 



a x <J v ~ h x a x — hy(Ty 



(2.4) 



Since we are studying isospin symmetric matter, we 
do not distinguish the up and down quark sectors. The 
mesonic sector has six parameters, that are the mass and 
couplings m 2 , Ai, A2 and c and the explicit chiral sym- 
metry breaking terms, h x and h y . They are adjusted to 
the pion and kaon decay constants f v and fx and meson 
masses of the scalar and pseudoscalar octet, m n , ttik-, 
m^ + rn^, and m a . The mass of the sigma meson is still a 
poorly known number, but the most recent compilation 
of the Particle Data Group [46| considers that m a can 
vary between 400 MeV and 550 MeV. Once given this 
set of masses and decay constants, the model parameters 
are defined. For their explicit expressions, please refer to 
Refs. [H, |4l[ . The values of the constants we use to cal- 
culate these parameters are listed in Table |U Wc will be 
using frequently the value m a = 500 MeV for the sigma 
mass. 

The second term in Eq. (|2.3[) is responsible for includ- 
ing the physics of color confinement through the intro- 
duction of a potential energy for the expectation value 
of the Polyakov loop. The explicit functional form of the 
potential energy density for the Polyakov loop is still not 
known directly from first-principle calculations. Instead, 
a common approach is to choose a functional form for the 
potential that reproduces crucial features of pure gauge 
theory and then adjust a set of free parameters to the 
results for the Polyakov loop and the thermodynamical 
observables of Monte-Carlo lattice calculations. A possi- 
ble parametrization of the Polyakov loop potential is the 
polynomial parametrization 0, [l~8| , 
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poly 



b2 (T) 



rp4 



G 
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(2.5) 
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TABLE II. Parameters of the gauge potential parametriza- 
tions for fits to the lattice simulation [471 ]. 



poly 


a 


ai 


a 2 


«3 


63 


64 


[12 


1.53 


0.96 


-2.3 


-2.85 


13.34 


14.88 


[18] 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 


log 


A 


A 1 


A 2 


B 3 






(19] 


3.51 


-2.47 


15.2 


-1.75 







with the temperature-dependent coefficient 62 defined as 

fT \ 3 

■ a 2 I — ) + a 3 



62 (T) = a + a 1 



(2.6) 



An equivalent polynomial parametrization had been pre- 
viously proposed in fl7| , with different definitions for the 
coefficients in Eq. (|2.5[) . A simple calculation allows the 
translation from one set of coefficients to the other. As 
far as the polynomial parametrization is concerned, we 
shall stick to the form (|2.5I) . 

Another possible parametrization for the effective po- 
tential of the Polyakov loop is provided by (i~9j 



It,, 



1 



J>4 

x In 



A (T) $$ + B (T) x 



1-6 +4 $ J + $' 5 -3 



(2.7) 



where both coefficients are temperature dependent 

2 



( y) +A2 (?) 



( Y 



(2.8b) 



The form (|2.7I) of the Polyakov loop potential is called 
the logarithmic parametrization. In Refs [r^ - [l9| the pa- 
rameters of both the polynomial and of the logarithmic 
parametrizations (see next) were adjusted to the lattice 
simulation of Ref . [13] ■ They are listed in Table [TTJ 

Originally, the parameter To was devised to corre- 
spond to the transition temperature in pure Yang-Mills 
theory To = 270 McV. However, in full dynamical 
QCD, fcrmionic contributions and the matter backreac- 
tion modify the pure gauge potential to an effective glue 
potential. Reference [6j estimated this running coupling 
of QCD by consistency with hard thermal loop pertur- 
bation theory calculations [H, |4t|. They mapped this 
effect to an ^/-dependent modification of the expansion 
coefficients of the Polyakov loop potential that results 
in a 7V/-dependence of To. The actual value of To for 
2+1 quark flavors with a current strange quark mass of 
95MeV [i| is T = 182 MeV. 

The authors in Ref. also estimated the dependence 
of the glue potential with the quark density and mapped 
it to a quark chemical potential dependence of Tq. Such 
a dependence can be expected in view of a //-dependent 
color screening effect due to quarks. The /1 dependence of 



To suggested in [6j is implemented as a small correction 
to the running coupling. Here, we will simply take a more 
simplistic approach, but without loss of generality, and 
consider To as a constant parameter. 

An improved mapping between the pure Yang-Mills 
potential and the quark-improved glue effective potential 
in full QCD is discussed and applied in Refs. [lO, [TTJ and 
will be taken into account in a future work. 

Finally, the last term of Eq. (|2.3[) represents the con- 
stituent quark sector coupled to the gauge field (repre- 
sented by the Polyakov loop variables) and to the mesons, 



n qi! (c- x ,a y ,$,$) = -2T E / 



d 3 p 

r 



f=u,d,s ' 

1 + 3 (<S> + §e- ( - E «'f~Vf )/T ) > 



x jln 1 + 3 ( 

xe -(E q ,f-^f)/T + e -3(E q , f -n f )/T~\ + 

+ In 1 + 3 ($ + <Pe- {E « ./+"/)/ T ) x 

Xe -(E q ,f+»f)/T + e -3{E qJ+IJ , f )/T I ^ 



(2.9) 



where (if is the quark chemical potential for each of the 
quark flavors. 

The constituent light and strange quark masses are, 
respectively, 



and 



(2.10) 



To fix the Yukawa coupling, we choose the constituent 
quark mass of the light (up and down) quarks to be mi = 
300 MeV in the vacuum, where (cr x ) = f n . This results 
in m s ~ 417 MeV for the constituent strange quarks. 

Notice that the dependence of fl q q with the chiral con- 
densates a x and a y in Eq. (|2.9I) is implicit in the quasi- 
particlc dispersion relation for the constituent quarks 



E qJ = Jk 2 



(2.11) 



One last comment should be made about the vac- 
uum (or sea) terms of the in-medium effective poten- 
tial. Many comparative studies within the NJL and lin- 
ear sigma models and comparisons with functional cal- 
culations [H, [50I l5lj have shown that vacuum terms can 
change qualitatively the structure of the T — fx phase 
diagram. In particular, if the sigma meson mass is suffi- 
ciently large, the line of first-order phase transition can 
eventually disappear due to the effect of the fermionic 
vacuum contribution at one-loop pol I52I [53| . In the 
present work, although we neglect the vacuum terms 
from both fermionic and bosonic fields, we adopt mainly a 
low value of the sigma mass, m„ = 500 MeV. Such a low 
value has been shown to allow for a first-order transition 
at low temperatures even with the inclusion of fermionic 
vacuum contributions only (see for instance Ref. Q and 
references therein). Therefore, we believe that the no-sea 
approximation does not spoil the general conclusions of 
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this work. Still, we agree with the authors of [53| that 
all the contributions from all fields of a model should be 
consistently taken into account, that is the case in func- 
tional calculations HO, HH . In this direction, a full study 
of the renormalized one-loop contributions from quarks 
and mesons in the PQM model will be taken into account 
in a future work 1541. 



C. Thermodynamical equilibrium 

Given the temperature T and the quark chemical po- 
tential /!, the effective potential (|2.3p is then given as a 
function of the four order-parameters of the model, a x , 
cr y , $ and $. In thermal equilibrium, the field configura- 
tions that contribute the most to the partition function 
are those that minimize the in-medium effective potential 
(or the in-medium effective action, for non-homogeneous 
systems). All other extrema of the effective action are ex- 
ponentially suppressed and give negligible contributions 
to the equilibrium thermodynamics of the system. 

Even though this issue has already been addressed be- 
fore in the literature (e.g. in Refs. [ID, EH)! we believe 
that it should be discussed more thoroughly, especially in 
the context of the PQM model. Our main motivations for 
pressing on this point are, first, theoretical consistency, 
second, that the results from the two approaches (the 
saddle-point approach and the present one) are different, 
and third, that only with equilibrium states described by 
minima of the effective potential can one calculate quasi- 
equilibrium properties of the system, such as the surface 
tension in a first-order phase transition, as we already 
commented above in Sec. [I] 

The one-loop effective potential of the PQM model at 
finite temperature and quark chemical potential, as de- 
fined in Eq. (|2.3[) . is a complex function of complex vari- 
ables. Therefore, it can have no minima and, as it stands, 
it cannot provide a standard description for a thermody- 
namical system in equilibrium. Finally, it makes sense 
to identify the extrema of a function with maxima, min- 
ima and saddle-points only if it is a real function of real 
variables. 

In order to make this statement clearer, let us explicitly 
identify the real and imaginary parts of the effective po- 
tential at finite T and Thus, we must write down the 
effective potential in terms of these real variables only. 

Let us start by making a change of variables in the 
potential by introducing the real and imaginary parts of 
the Polyakov loop variables as 

$ + $ $ — $ 
a = —— and (3 = . (2.12) 

2 2i v ' 



2 The following discussion will make it clear that at fj, = 0, $ and 
4> are identical and real, so that there is no sign problem in that 



The polynomial potential (12.51) and the logarithmic po- 
tential (|2.7j) , respectively, can be both rewritten as func- 
tions of the real variables a and f3 as 



and 

rp4 



-|(a 2 +/3 2 )-|(a 3 -3 a /3 2 )+^(a 2 +/3 2 ) 2 , 

(2.13) 



+B(T) log[l-6(a 2 +^ 2 ) + 8(a 3 - 3a/3 2 ) - 3(a 2 + f3 2 ) 2 ] . 

(2.14) 

We now turn to the contribution from the quarks at 
finite temperature and chemical potential. Once $ and 
$ are complex-valued variables, the potential (|2.9p can 
take complex values. Let us rewrite it as a function of 
real variables to make this statement explicit. 

Dropping the flavor indexes to keep the notation sim- 
pler, we define 

Z + EE 1 + 3($ + $ e -(^)/T )e -(B-M)/T + e -3(B-„)/T 

(2.15) 

and 

z_ ee 1 + 3(# + $ e -(^)/ T ) e -( J5 +^/ T + e-HE+rf/T^ 

(2.16) 

such that (12.91) can now be written as 



(2.17) 



After using (|2.12p and performing some straightfor- 
ward manipulations, one can see that the argument of 
the logarithm in (|2.17[) is complex, that is, 



= R + il , 

where 

R =1 + e -3(B-M)/T + e -3(E+n)/T + e -6E/T + 



(2.18) 



+6ae 



-E/T 



cosh 



+6ae~ iE / T 



(|)+ e - £/T cosh(^ 



2E/T\ e -2E/T 



+9(a 2 + /3 2 )(l + e 
+ 18(« 2 -/3 2 )e- 3 ^cosh(^) , 



(2.19) 



and 



I =6/3e- £ / T 



sinh 



T 



+6/3e 



-4E/T 



e- £/T sinh 



sinh (I 



(I) 



-36a/3sinh(^)e- 3B / T . 



(2.20) 
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The complex argument of the logarithm can be written 
in polar form, R + il = pe , with 



p = \JR 2 + I 2 and 



arctan | — 
K 



(2.21) 



so that the potential can be cast in a manifestly complex 
form, 



where 



"qq — "qq \ <■ J 'qq , 



(2.22) 



and 



^---2^E/^ ( 2 - 24 ) 



The imaginary part (|2.24[) is the manifestation of the 
fermion sign problem in the context of the PQM model 
already at the one-loop level. Note that it is very closely 
related to the sign problem in the PNJL model as well 
[2, 0, 0- An important aspect of Eq. (|2~2"4"|) is that it 
vanishes for p — 0, so that the effective potential becomes 
real and free of the sign problem. Furthermore, (|2.24p is 
odd in /?, while the real part (|2.23p is even in (3. This 
means that we must have (f3) = 0, i.e., ($) = ($) for 
p = 0, which is a well-known result. 



D. A Comment about the Sign Problem 

In order to make the sign problem explicit in the PQM 
model, let us write the grand partition function for the 
model in the mean-field approximation, 



Z 



D*] cxp 



(2.25) 



where V is the volume of the space and ^ = (a x , a v , a, j3) 
formally represents the (real) order parameters. 

Notice from Eqs. (|2~20f and (|2~24|) that the imaginary 
part of the quark contribution is odd in (3 and all the 
other contributions arc even in j3. Let us then split the 
effective potential in a /3-odd part (that coincides with 
Qqq) and a /3-even part. The functional integral is to 
be performed for every possible (real) value that ft can 
assume. We can organize the sum such that the contribu- 
tions for a given (3q and its negative — /?o are assembled. 
After a few simple manipulations, one finds 



Z= I [D^] exp 

hem. 



— (ftp-even + ^qq) 



[D*] exp 



" V 






"^^(3 — even 


COS 


y"qq 



(2.26) 



The integrand in the expression (|2.26j) for the parti- 
tion function is not positive defined, as a sound Boltz- 
mann factor should be. Recall that the integrand of the 
partition function (the density matrix elements) corre- 
sponds to a sort of probability density, which must be 
non-negative. This is not true for the grand-partition 
function (|2.26|) and we conclude that the PQM model 
has the sign problem at the mean-field level for a finite 
chemical potential. 

If we insist on writing the partition function (|2.26|) in 
the same form as (|2.25p . we end up with 



T 

S2 = £lp- even — — log 



V r 
I y^qq 



(2.27) 



where V is the volume of space. This function is not 
physically acceptable as an effective potential. First, it 
is a volume-dependent effective potential (therefore, it is 
not intensive) and second, it is not defined in the ther- 
modynamical limit V — > oo. 

One possibility to circumvent the sign problem, as in- 
dicated in [l5| is to treat the imaginary part of the effec- 
tive potential pcrturbatively in an expansion in powers 
of T/V. In the first order of the approximation (which 
the authors of [l5[ identify with the mean-field approxi- 
mation), one simply ignores the imaginary part |2.24p of 
the effective potential at finite chemical potential How- 
ever, if one neglects f2^, the expectation value of f3 is 
zero due to the even parity of with respect to j3. As 
a consequence, the difference (<& — <&) = 0, which is in dis- 
agreement with complex Langevin [55j and Monte-Carlo 
simulations [J . In spite of this setback, we understand 
that our approach is simply an approximation scheme 
that has the strong theoretical advantage of dealing with 
the well-established minimization procedure for finding 
the state of equilibrium. 



E. Finding the Minima of the Effective Potential 

From now on, we shall neglect the imaginary part of 
the effective potential, as explained before, i.e., consider 

n = U{a x ,a y )+U(a,0) + n q \(a x ,a y ,a,p) . (2.28) 

Given the temperature T and the quark chemical po- 
tential p, the effective potential (|2.28p is then given as a 
function of the four order-parameters of the model, cr x , 
(j y , a and /3. In equilibrium, the expected values for the 
order parameters are given by minimizing the effective 
potential (|2.28p with respect to each of them. There- 
fore, the first necessary condition is to find a point in the 
order-parameter space such that 

dn on dCl dCl , 

^ = ^ = ^r = ^ = ' 2 - 29 

oa x 0(7 y oa op 



3 In this model, ignoring is equivalent to taking the modulus 
of the Dirac determinant that is tacitly present in 112.261 1 
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when the derivatives are evaluated at the extrema. 

The second condition is the positivity of all the eigen- 
values of the 4x4 Hessian matrix 



as 



(2.30) 



with (X{ = o~ x , a y , a, 0) evaluated at the points where 
([2~29]1 is satisfied. 

Hence, we must calculate the first and second deriva- 
tives of the effective potential Ct with respect to the order 
parameters in order to find its minima. 

Moreover, once we are interested in phase transitions, 
our algorithm to find the minima considers the possibility 
of multiple minima. It performs the following basic steps. 
With a given starting position in the four-dimensional 
space of order-parameters, the algorithm looks for a point 
with vanishing gradient (|2.29l) . As a second step, the 
eigenvalues of the Hessian matrix (|2.30|) are evaluated at 
this point and their signs are checked: if they are all pos- 
itive, this point is a minimum and it is saved as such; 
otherwise it is either a saddle-point or a maximum and, 
therefore, it is discarded. Next, another starting position 
in the order-parameter space is tried, and the procedure 
is repeated. If the point found is a minimum, it is com- 
pared to the previous points: if all its four coordinates 
are sufficiently close to any of the other points (within 
1%), it is discarded; otherwise, it is saved as a new min- 
imum. By sweeping a sufficiently wide window in the 
parameter space, this is a simple algorithm able to find 
all the minima of the effective potential. 

The thermodynamical state of equilibrium for each 
given T and fi is given by the values of the order pa- 
rameters at the global minimum of the effective poten- 
tial, found after the prescription described above. This 
allows us to describe any thermodynamical quantity of 
equilibrium. Particularly, we can study the behavior of 
the order parameters as functions of the temperature and 
chemical potential, as well as the phase diagram of the 
model. Moreover, finding the minima of the effective 
potential at the vicinity of a first-order phase transition 
allows the description of the dynamics of the transition. 
More specifically, one can calculate the surface tension 
between the two phases, as we discuss in the next sec- 
tion. Our results will be discussed in Sec. IIVI 



F. The Saddle-Point Approach 

In the literature dedicated to the thermodynamics of 
strongly interacting systems described by chiral models 
coupled to the Polyakov loop, such as the PNJL and 
PQM models, the most common method to determine 
the state of thermodynamical equilibrium is sometimes 
called the saddle-point approach. It consists of finding 
the points in which the in-medium effective potential of 
the theory has vanishing derivatives with respect to the 
order parameters, as in Eq. (|2.29l) . These extrema, how- 
ever, are not usually required to be minima of the ef- 



fective potential. Actually, the effective potential often 
has no minimum whatsoever as a consequence of one of 
the two following reasons. Either the effective potential 
is considered in full, so that it is a complex function of 
complex variables (as we have discussed previously), or 
the Polyakov loop variables arc restricted to being real 
(instead of complex) numbers (even before their expecta- 
tion values are evaluated). In the first case, no minimum 
can be defined because the effective potential is complex 
valued. In the second case, it can be easily seen that the 
effective potential is unbounded from below^ @. 

The approach we propose in this section leads to pre- 
dictions that differ from those made by the saddle-point 
approach at \i ^ 0. Some of these differences will be 
briefly discussed later on in Sec. IIVI 



III. NUCLEATION IN THE PQM MODEL 

A. Homogeneous Thermal Nucleation 

The dynamics of a first-order phase transition at 
small metastability can be described by phenomenologi- 
cal droplet nucleation models [2(], [56l - l59| . In such a fam- 
ily of models, the transition between a metastable and a 
stable phase takes place by the appearance and growth of 
domains (droplets or bubbles) of the stable phase inside 
the metastable phase. The phase conversion is finished 
when these domains grow and coalesce completely. In any 
case, a minimum-sized bubble is needed for the beginning 
of the phase transition, as can be inferred from the fol- 
lowing heuristic argument. The bulk free energy density 
of the metastable phase (often called false vacuum) is, by 
definition, higher than that of the stable phase (the true 
vacuum). Therefore, the conversion of a given fraction 
of the system into the stable phase makes the bulk free 
energy of the whole system lower. However, given that 
such a conversion takes place within a connected domain 
of the system (most likely in a spherical bubble [poj). 
an interface is needed in order to separate the (stable) 
interior from the (metastable) exterior of this domain. 
Once the creation of an interface represents an energy 
cost, the mechanism of phase conversion through bubble 
nucleation settles a competition between the free energy 
gain from the phase conversion of the bulk and the energy 
cost from the creation of an interface. Roughly, one can 
say that the free energy shift due to the appearance of a 
spherical bubble of the stable phase of radius R inside a 
metastable system is 



AK = 



4~ 
3 



-^-R 3, /stable + 47Ti? 2 / waU 



4~ 
3 



R /metastable 



4 Just consider $ = in Eq. I I2.5H or in Eq. 112. 7t and we can realize 
that the Polyakov loop potential is unbounded from below for 
§ — > oo. 
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y# 3 A/ + 47rE 2 /wall, 



(3-1) 



where /stable and /mctastabio are, respectively, the bulk 
free energy densities of the stable and metastable phases, 
and / W aii is the surface energy density of the bubble wall, 
that is, the surface tension of the interface between the 
two phases. This formula clearly shows the competition 
between bulk (negative) and surface (positive) contri- 
butions. Notice that the shift in the bulk free energy 
(A/ = /stable - /metastable < 0) is proportional to the 
volume of the bubble, while the surface free energy cost 
is proportional to its area. For the nucleation of small 
bubbles, the energy cost is higher than the energy gain. 
Therefore, small bubbles shrink. On the other hand, a 
very large bubble represents a large bulk energy gain, 
which is higher than the surface energy cost in creating 
the bubble. As a consequence, large bubbles tend to grow 
even more and to occupy the whole system, completing 
the phase transition. Consequently, this energy compe- 
tition implies the existence of a so-called critical bubble: 
any bubble smaller than the critical bubble will shrink 
and any larger bubble will grow and drive the phase con- 
version. For this reason, the critical bubble is the crucial 
object in the theory of dynamical first-order phase tran- 
sitions of slightly metastable systems. 

The appearance of a bubble (critical or not) of the sta- 
ble phase inside a metastable system is a natural conse- 
quence of the never-ending thermal and quantum fluctua- 
tions of any thermodynamical system sufficiently close to 
a first-order phase transition. As we just discussed, each 
bubble created by these fluctuations may grow or shrink, 
depending on its energy budget with regard to a homo- 
geneous metastable phase. One should also have in mind 
that larger fluctuations (like a critical bubble) should be 
less common than smaller ones. Although small bubbles 
are frequently created, they rapidly disappear and do not 
contribute to the process of phase conversion (with the 
exception in a weak first-order phase transition, when co- 
alescing subcritical fluctuations [HJ HH can complete the 
phase transition without the nucleation of critical bub- 
bles) . Only those that have a size equal to or larger than 
the critical bubble have a decisive role. The smallest (and 
therefore the most probable) among them is the critical 
bubble. This means that the mean time that it takes 
for random fluctuations to create a critical bubble is the 
shortest time scale for the creation of a lasting domain 
of the stable phase, which is the dynamical seed of the 
phase conversion. 

Let us assume that the system is in a metastable state 
in quasi-equilibrium with a reservoir with intensive coor- 
dinates generically represented by 3? (e.g., temperature, 
chemical potential, etc). Being metastable, bubbles of 
the stable phase with different sizes randomly appear and 
subsequently disappear. This process keeps happening 
until a critical bubble is nucleated and the phase con- 
version effectively starts. It can be shown by different 
approaches that the rate at which critical bubbles are nu- 
cleated per unit time, per unit volume can be expressed 



in the form [20|, |21|, |56|, |63|-[6J 

T(X) = y(T) exp 



AF&(3£) 



T 



(3.2) 



where T is the temperature of the system in equilibrium 
with the reservoir. The pre-exponential factor (or pref- 
actor) J'(IR) corresponds to the probability for a critical 
bubble-like field fluctuation </>b to be generated and grow 
[20l - [22j . The last factor in (|3.2I) is a Boltzmann factor in 
which AFh (31) is the shift of free energy (as compared to 
the homogeneous metastable phase) due to the formation 
of a critical bubble. It can be easily shown that AF{,(3?) 
can be cast as in (|3.ip for a small degree of metastabil- 
ity, where the thin-wall approximation is valid [6r| and 
where it is proportional to /^ all /(A/) 2 . In spite of the 
general importance of the prcfactor in (|3.2[) . its specific 
form is not crucial for the nucleation rate at a small de- 
gree of metastability. Close to the coexistence of both 
phases, it shows a proportionality to / wall /A/, so that 
the nucleation rate is strongly dominated by the expo- 
nential factor [2l|. For this reason, it will be enough to 
focus on the free energy shift in this work. 

The process of bubble nucleation in an impurity-free 
environment is called homogeneous nucleation. In this 
work we shall only consider the process of homogeneous 
nucleation, which is not the most common in natural en- 
vironments, such as in a boiling liquid. In such cases, 
the presence of impurities can drastically accelerate the 
nucleation of bubbles, and the process is called inhomoge- 
neous nucleation (which is also the case when subcritical 
thermal fluctuations can dominate [6l|, H3)- The process 
of inhomogeneous nucleation can be orders of magnitude 
faster than homogeneous nucleation because impurities 
(like dust) often reduce the free energy cost for the for- 
mation of a critical bubble, raising the probability for its 
formation [6?], [68[ . We do not consider it in this work for 
two reasons. First, we wish to underestimate the nucle- 
ation rate and inhomogeneities could only increase this 
rate, and second, we wish to keep the approach as simple 
as possible. 



B. The Coarse-Grained Free Energy for a Single 
Scalar Order-Parameter 

As discussed in 0, [H, [H [70[ , the nucleation rate of 
critical bubbles can be calculated from the microphysics 
using scmiclassical methods in Euclidean thermal field 
theory. We stress the importance of considering the ef- 
fective action in the problem of bubble nucleation, and 
not simply the effective potential, once a critical bub- 
ble is clearly a non-homogeneous field configuration. For 
simplicity, let us start with the Euclidean Lagrangian 
density for a single scalar order-parameter field (which 
we generically call (f>) of the form, 



.(d^f + V{<j>). 



(3.3) 
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In this simple model, one assumes that the order parame- 
ter for a system in thermodynamical equilibrium is given 
by the expectation value of <f). In general, it depends on 
the properties of the reservoir, such as its temperature or 
chemical potential, which we generally denote by 31. The 
(Euclidean) action is 



dr \ d 3 x£ s [0(x, r)]. 



(3.4) 



In the high-temperature limit, j3 = 1/T — > 0, the imag- 
inary time dependence of the order parameter can be ne- 
glected [(39[ and, therefore, we make the approximation 



S E [<j>,K] 
where one identifies 



T ' 



d 6 x 



(V^) 2 + V eS (c/>, 31) 



(3.5) 



(3.6) 



with the coarse-grained free energy of the system. Notice 
that the coarse-grained free energy is a sort of dimension- 
ally reduced effective action, so that the tree-level poten- 
tial V{4>) must be replaced by the medium-dependent 
effective potential V^ff (0, 31). In full thermodynamical 
equilibrium, the minimization of the coarse-grained free 
energy (which is equivalent to the minimization of the 
Euclidean action) is achieved by a constant field configu- 
ration 0(x) = (/)q so that the gradient term vanishes and 
V e g((f>Q, 31) must be a global minimum of V e ff. 

The possibility of a mctastable state arises when V e g 
develops some local minimum other than the global min- 
imum at (j) = <^o- In this framework, a metastable state 
is described by a constant field configuration 0/ that is 
a local minimum of V e g. For this reason, this second 
minimum is often called a false vacuum of the potential, 
while the global minimum is called the true vacuum of 
the theory, (f> t = 4>o- A bubble is then represented as a 
non-homogeneous spherically symmetric field configura- 
tion <f)(r) such that [H, H3, Q 



lim 



p ) = 4>f and 



(3.7) 



where (f>f is the value of the order-parameter field at the 
false vacuum. That is, away from the center of the bub- 
ble, the system is in the metastable phase. But, in the 
vicinity of its center, the field configuration should be 
close to the stable minimum (but not necessarily exactly 
on it). 

The critical bubble is a saddle point field configuration 
4>b that extremizes the functional F, i.e., it solves the 
Euler-Lagrange equation 



8F\6- T] 







V 2 0(x) 



dV cf 
d<t> 



[0(x)]=O. (3.8) 



It can be shown [60( for a wide class of Lagrangians, 
including (|3.3[) , that the smallest value of F indeed corre- 
sponds to a spherically symmetric solution of (|3.8p . 4>(r), 



so that the equation to be solved is now the non-linear 
ordinary differential equation 



(r) , 2 d<t>(r) _ dV cH 



dr 2 



r dr 



Mr)), 



(3.9) 



with the boundary conditions (|3.7I) . 

The coarse-grained free energy associated with a spher- 
ical bubble (j)b is then 



Fh = 4tt 





d(f>b(r) 




dr 



+ V eS [Mr)} } , (3.10) 



which directly follows from Q3.6p . 

Once the solution for (|3.9[) and (|3 . T[) is found, that is, 
the critical bubble profile 4> C rit(r), the shift in the coarse- 
grained free energy due to the appearance of a critical 
bubble, 



AF b (3l) = F[ 



31] 



F[<j> f - 31], (3.11) 
, can be readily calcu- 



needed in the nucleation rate 
lated. 

For a generic effective potential V e g, the solution of 
(|3.9p with boundary conditions (|3.7[) cannot be obtained 
analytically. However, an approximate solution can be 
found when the system is very close to the coexistence 
line, so that it is slightly mctastable and the thin-wall 
approximation [63|, [6g, E^] is applicable. Within these 
limits, the coarse-grained free energy shift p. lip can be 
well approximated by the expression 



AF b (3l) 



16tt S 3 
~(AV4ff) 2 ' 



where 



S(T) 



/ dr 


"#critW 


'o 


dr 



(3.12) 



(3.13) 



is the surface tension of the critical bubble interface be- 
tween the phases. Notice that the surface tension is cal- 
culated directly from the critical bubble solution <fi cr it. 
It must be so because the surface tension must con- 
tain information about how the system reacts to inho- 
mogeneities (e.g., a wall). That is, any description of a 
critical bubble has to take into account more than just 
the bulk thermodynamics. This is the reason why the 
coarse-grained free energy (|3.6p is needed in the formal- 
ism for bubble nucleation. 

The quantity AV e f{ = V^^t) — (</>/) is the differ- 
ence between the bulk free energy in the two homoge- 
neous vacua. In the grand canonical potential, it can 
be identified as AV c ff(T,^i) = — Ap(T, /i), i.e., minus the 
difference of pressures between the two phases. Using 
the thin-wall approximation, the surface tension integral 
(|3.13p can be calculated without solving for the profile. 
After changing variables from r to </>, one finds 



r <j> 

W) = \ 

J<t>t 



(3.14) 
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so that only the effective potential 14ff[</>; Ut] is needed 
to calculate S(!R) in the thin- wall approximation. Notice 
that V e g is normalized so that its global minimum is lo- 
cated at 4> = 4>t = 0- Exactly at the coexistence points, 
where the thin-wall approximation is exact, the minima 
are degenerate with V e s = 0. 

As a last remark, let us notice that the surface tension 
cannot be correctly defined unless <pt and (f>f are actual 
minima of the effective potential. This a further moti- 
vation for the careful approach discussed in the previous 
Sec. El 



C. The Coarse-Grained Free Energy for the 
Polyakov-quark-meson Model 

In order to calculate the free energy shift AFt(3Vj due 
to the nucleation of a critical bubble within the PQM 
model, we first need to define the coarse-grained free en- 
ergy functional and then identify an order parameter. 

Just as in the case of a single order-parameter, the 
coarse-grained free energy of the PQM model has its ori- 
gin in the in-medium effective action of the theory. In the 
PQM model, the Lagrangian density of the chiral fields 
<j x and <j v directly leads to a kinetic term of the form 
|V^>| 2 . The kinetic term for the Polyakov loop variable, 
however, is not determined a priori. For simplicity, once 
we consider the Polyakov loop order parameters a and /3 
as independent, real variables at finite /i, we assume the 
kinetic term 



£kinK/3] 



^(M 2 + ^(d,/3) 2 , 



(3.15) 



which is motivated by the Z(3)-line model of Ref. j7l| . 

Effective models cannot do much more than estimate 
the value of the kinetic parameter re. From dimen- 
sional arguments, it is estimated to be re 2 = NcT^/gl 



[171 |71| . Assuming N c = 3 and 



2/4tt ~ 0.3, 



we find re ~ 0.9To, which is of the order of magnitude 
of the only scale in the pure glue model, the transi- 
tion temperature Tq. Another consistent approach is 
to consider as an input the surface tension of the pure 
gauge SU (3) theory calculated through lattice Monte- 
Carlo simulations [721 [73j and in an effective matrix model 
[74[, Ss(7(3) — 0.016Tq. With a given parametrization of 
the Polyakov loop potential, the parameter re can be fit- 
ted from the value of the surface tension at T = Tq using 
Eq. (|3.14[) . We adopt this second approach here in order 
to fix the parameter re with the implications discussed in 
Sec. El 

The coarse-grained free energy for the PQM model can 
be written as 



where £l(cr x , <J y , a, /3) is the effective potential defined in 
Eq. (l2~28l) . 

As discussed previously, the critical bubble is a solution 
of the four coupled Euler-Lagrange equations that arise 
from the functional (|3.16[) . There is no general procedure 
to solve this set of coupled equations, even though some 
ansatz solutions can be eventually tried for very simple 
effective potentials [75j . 

There are two possible ways in which we can tackle the 
problem. The first is, of course, to numerically solve the 
four equations that follow from the extremization of £F 
simultaneously. The exact solution will define a path in 
the four-dimensional space of order parameters. Notice, 
however, that this path is in general not in the "valley" 
that connects the two minima, as can be intuitively seen 
from the inverted potential mechanical analog. 

Let us now argue that our approach will lead us to 
an overestimate of the surface tension or, equivalently, to 
an underestimate of the nucleation rate. According to 
Eq. (|3.2[) . the higher the free energy shift of the critical 
bubble, the lower is the nucleation rate (at least as long 
as the thin- wall approximation is valid). Now consider 
that the exact critical bubble solution is found. Consider 
also a small deviation from it that follows a different path 
in the order-parameter space. Once the true critical bub- 
ble is a saddle point solution in a functional space, one 
expects that the distorted path will have a higher value 
of AF as compared with the true bubble. Hence, we can 
(artificially) constrain the configuration path to a given 
arbitrary line that connects the two vacua in the space of 
order parameters. This path will give us an overestimate 
of the free energy and, therefore, of the surface tension. 
The simplest choice is, of course, a straight line that con- 
nects both minimal! For example, let a x and a x 2 ^ be 
the values of the <j x order-parameter in the two minima 
of Cl close to the coexistence line. The interpolation 

a x = + (1 - £)4 2) (3.17) 

is such that for < £ < 1, the value of a x varies from 
one minimum to the other. If the same function £ is used 
to parametrize the path followed by the remaining order 
parameters, 

= + (i-o4 2) < 

a=5« (1) + (l-()a (2) , 

/3 = £/3W + (l-£)/3 (2) , (3.18) 

then this path is a straight line in the four-dimensional 
order parameter space. 

It is now natural to define the four-dimensional order- 
parameter for the PQM model as 



* = {cr x ,a y , rea, re/3) . 



(3.19) 



J[a x ,a yi a,p\ 

.2 



d 3 x 

2 



^(Vcr^ + ^V^) 2 - 



, (3.16) 



5 In the whole range of interest in the T — fj, plane, we have always 
found either one or two minima. The case of three or more 
minima, which we have not encountered, would lead to a much 
more complicated analysis. 
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It can be easily shown from Eqs. p,17p and (|3 . 19[) that 
if we write the coarse-grained free energy (|3.16p in terms 
of the field £, it assumes the form 



1 ^ 



HO = / d A x 



Vh2 



(vo 2 +n(o 



(3.20) 



where 



h z = (* - ^ 

= (Aa x ) 2 + (Aa y ) 2 + («Aa) 2 + (kA/3) 2 , (3.21) 

and fi(£) is the projection of the effective potential O(^) 
along the straight line defined by (|3.17[) . Notice that 
the coarse-grained free energy 2r is formally equivalent 
to the single field coarse-grained free energy (|3.6[) and 
we can now consider £ as a scalar order-parameter. As 
a result, the surface tension in the PQM model can be 
overestimated by 



(3.22) 



where the domain of integration ranges from one mini- 
mum of (and thus of f2) to the other along a straight 
line path in the four-dimensional order-parameter space. 

Notice that a path other than a straight line in the 
space of order parameters would not lead to a simple 
coarse-grained free energy with the same form as (|3.20[) . 
In the general case, even the definition of the surface 
tension should be reformulated. This problem will be 
treated in a separate work. 

In the following section, we present our results for the 
equilibrium thermodynamics as well as for the surface 
tension in the PQM model. 



IV. RESULTS AND DISCUSSION 

We will now present our original results for the Nf = 
2 + 1 PQM model for quantities in strict thermodynami- 
cal equilibrium, as well as an overestimate of the surface 
tension of a hadron-quark interface as described by this 
model. We use the methods and approximations dis- 
cussed in the previous sections. Here, we call the Poly- 
nomial parametrization of [13] "Poly-I", while that of 
Ref. [18l| is called "Poly-II" . The logarithmic parametriza- 
tion of [19[ is referred to as "Log" . 

Notice that the imaginary part of the effective poten- 
tial (|2.24j) is zero for \i = 0, and our scheme has to 
be completely equivalent to the saddle-point approach in 
this limit (notice that a = $ for /i = 0). Our numerical 
results show that this is true indeed. 

In an explicit crosscheck, we calculated the pseudo- 
critical temperatures for deconfinement and chiral sym- 
metry restoration at /x = 0, with the same parameter 
set of Ref. 0. The chiral pseudocritical temperature 
is defined as the peak in the chiral susceptibility \x = 



0.8 
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FIG. 1. (Color online.) Normalized nonstrange condensate 
{{<T x )/fn, solid line) and strange condensate {{cr y ) / f-n, dashed 
line) and Polyakov loop variable ({a), dotted line) as a func- 
tion of the temperature for /i = 0. The parameter set used 
is that of Table Q] with m CT = 600 MeV and T = 270 MeV 
together with the logarithmic Polyakov loop potential. 



—d(cr x }/dT, while the deconfinement pseudocritical tem- 
perature takes place at the peak of \p = d(a)/dT. We 
found the same pseudocritical temperatures as those in 
for the Poly-I and the logarithmic parametrizations 
at /i = 0. This shows numerically a consistency between 
our approach and the saddle-point approximation in the 
regime where both methods have to lead to the same re- 
sultsE However, we must stress, the two approaches do 
give different predictions at finite quark chemical poten- 
tial, as expected. This difference is mainly due to the 
fact that we are explicitly considering only the real part 
of the effective potential at finite chemical potential. As 
we have already discussed in Sees. UlandHIl this is a nec- 
essary condition in order to fulfill the need of defining 
minima and for their existence in the thermodynamical 
potential, which can only be achieved through the real- 
ity of the effective potential. In the following, we present 
our original results alongside the corresponding calcula- 
tions we performed using the saddle-point approach for 
comparison. 

Let us first discuss our results for the behavior of the 
order parameters for the model at fi = for the loga- 
rithmic parametrization of the Polyakov loop potential. 
The behavior for the polynomial parametrizations Poly-I 
and Poly-II are very similar and we do not display them 
here. Our choice of parameters is shown in Table HI with 
m a = 600 MeV and T Q = 270 MeV. In Fig. [TJ we show 
our results for the expectation value of the chiral con- 
densates o~ x and o y (normalized by the vacuum value 
of the nonstrange condensate {p x ) vac = f w ), as well as 



6 We refer the interested reader to Ref. [tJ for the actual values of 
pseudo-critical temperature found within the Nf =2 + 1 PQM 
model. 
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FIG. 2. (Color online.) Normalized nonstrange condensate 
{{(fx)/ fn, solid line) and strange condensate ((fy)//^, dashed 
line) and Polyakov loop variable ({a), dotted line) as a func- 
tion of the chemical potential for T = 2 MeV. The param- 
eter set used is that of Table [J with m a = 500 MeV and 
Tq — 182 MeV. For the Polyakov loop potential the logarith- 
mic parametrization is chosen but at this low temperature its 
influence is negligible. 



for the Polyakov loop variable a. The steep change of 
the order parameters at around T = 205 MeV indicates 
the pseudo-critical temperature of both the chiral and 
the deconfinement crossovers. (For some other choices of 
parameters, these pseudo-critical temperatures may not 
coincide Q.) 

The low-temperature behavior of (o~ x ) and of (a y ) 
shows little sensitivity to the choice of the Polyakov loop 
potential, being dominated by the chemical potential- 
dependent contribution from the mesonic sector. 

Next, we observe that, at low temperatures, the or- 
der parameters have a discontinuity at some finite /z. 
This signals a first-order phase transition. In Fig. [5] we 
show the behavior of the condensates <r x and a y and the 
Polyakov loop variable a (multiplied by a factor of ten) as 
a function of the chemical potential for a very low tem- 
perature T = 2 MeV with sigma mass m CT = 500 MeV. 
Notice that the Polyakov loop variable has a very low 
expectation value, which grows very slowly at such a low 
temperature, even at high densities. 

Having found a crossover at \i = and a first-order 
transition at T = 0, the next step is to calculate the 
complete T — ji phase diagram of the model. The phase 
diagram is shown in Fig. [3] for m a = 500 MeV and To = 
182 MeV for the three parametrizations of the Polyakov 
loop potential we have considered in this work. 

As mentioned before, the saddle-point approach leads 
to different results to quantities in thermodynamical 
equilibrium. Comparing the phase diagrams obtained 
within the two methods, we notice that they are equiv- 
alent at \x = (where the sign problem does not exist) 
and at T = (where only the chiral fields are relevant 
and, again, the sign problem does not exist). For this 




Quark Chemical Potential [MeV] 



FIG. 3. (Color online.) Lines of first-order phase transition 
of the PQM model. The sigma meson mass is m CT = 500 MeV 
and To = 182 MeV. Parametrizations of the Polyakov loop 
potential: Poly-I (full black line), Poly-II (dashed red line) 
and Log (dash-dotted blue line). 



TABLE III. Temperature and quark chemical potential co- 
ordinates of the critical end point (CEP) found with the 
minimization procedure of Sec. [Ill (Tcfc), an d with the 
saddle-point approach, (Tc,^c) S - We choose a sigma meson 
mass m a = 500 MeV and pure glue transition temperatures 
To = 182 MeV (with effective screening) or T = 270 MeV 
(pure gauge). The parameters not shown in this table are 
found in Table |U 



Parametrization 


To [MeV] 


(T c ,Mc)[MeV] 


(T c ,/xc) S [MeV] 


Log [19] 


182 


(143, 129) 


(143, 128) 




270 


(192,88) 


(192,84) 


Poly-I [17] 


182 


(139,99) 


(140,92) 




270 


(171, 103) 


(175,83) 


Poly-II [18] 


182 


(146, 115) 


(152,80) 




270 


(176, 129) 


(184, 103) 



reason, the phase diagrams can only differ between these 
two extrema and, in particular, they can lead to differ- 
ent positions for the critical end point (CEP). Table ILTTl 
shows the T — fx coordinates of the CEP found using the 
method we propose and the saddle-point method. De- 
pending on the parametrization and value of Tq the po- 
sitions of the CEP found within the two methods can 
differ significantly. The difference increases with a larger 
value of the glue transition temperature To and is most 
pronounced if using the Poly-II parametrization. Due to 
the slope of the phase transition line around the CEP 
the differences of the critical chemical potential is larger 
than that of the temperatures. 

The effective potential has always two degenerate min- 
ima over the coexistence line of the phase diagram, each 
one with a different quartet of expected order parame- 
ters. In Fig. |H we show how the coordinates of these 
minima evolve along the first-order line of the phase di- 
agram, starting from T = and finite /i up to the CEP. 
Notice that the two minima smoothly merge at the CEP 
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FIG. 4. (Color online.) Degenerate values of the chiral con- 
densates a x and a y (upper part) and Polyakov loop variable 
a (lower part) over the first-order line in the T — /j, plane as 
a function of the coexistence temperature. Notice that the 
minima merge smoothly at the critical end point. The pa- 
rameter set used is that of Table U with m CT = 500 MeV and 
To — 182 MeV together with the logarithmic parametrization 
of the Polyakov loop potential. 



with diverging derivative with respect to the tempera- 
ture. This signals the expected second-order phase tran- 
sition at the CEP. 

The two minima of the potential that are degener- 
ate on the coexistence line persist as a global and a 
metastable local minimum in some region of the phase 
diagram around the phase transition line. Going away 
from the coexistence line, the intervening maximum ap- 
proaches the local minimum until these two extrema meet 
and form an inflection point that defines the spinodal 
line. In Fig. [S] we show the extension of the metastable 
region of the T — fj, phase diagram that is limited by the 
spinodal lines. The extension is relatively independent of 
the parametrization of the Polyakov loop potential. We 
find that the degree of metastability that can be reached 
is relatively modest. 

Let us now study the behavior of the surface tension in 
the PQM model along the coexistence line for different 
paramctrizations of the Polyakov loop potential IX. The 
surface tension was evaluated in the thin- wall approxima- 
tion, according to Eq. (|3.14p . In Fig. HI we show the be- 
havior of the surface tension along the first-order transi- 
tion line for the Poly-I, Poly-II and the Log parametriza- 
tions of the Polyakov loop potential. The surface tension 
for the Poly-II potential is off scale to the other two ones 
and, therefore, it is shown in the inset of Fig. [5] 

Even though equilibrium observables like the behavior 
of the order parameters and the phase diagram are quite 
similar for the three parametrizations, there is a remark- 
able difference between them regarding the surface ten- 
sion, in particular for the result obtained for the Poly-II 
parametrization. We interpret this as an artifact of the 
Poly-II parametrization. Although all paramctrizations 
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FIG. 5. (Color online.) Metastable regions of the PQM 
model for the three considered parametrizations. We compare 
their extension with respect to the coexistence line (/x C oex)- 
The temperature is given relative to the one of the criti- 
cal end point Tc, see Table IIHI The sigma meson mass is 
m a = 500 MeV and T = 182 MeV. 
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FIG. 6. (Color online.) Surface tension relative to its zero- 
temperature value along the first-order transition line in the 
T — fj, phase diagram as a function of the temperature relative 
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Log parametrizations with To = 182MeV. Inset: The Poly- 
II parametrization is shown for comparison. The parameter 
set used is that of Table U with m CT = 500 MeV. The other 
parameters can be found in Tables [TTTl and ITVl 



lead to very similar descriptions of the thermodynam- 
ics of the dcconfincmcnt transition in equilibrium, they 
offer different behaviors for the order parameters. At 
the pure glue transition temperature To, the Polyakov 
loop potentials present two minima, which correspond to 
the two coexisting values of $ = $ at T = Tq. One of 
these minima is always at $ = 0, representing the con- 
fined phase, while the other sits at some nonzero value 
$ = $0- Using the parameter sets shown in Table ILTl 
one can easily see that the Poly-I parametrization pos- 
sesses the second minimum at $0 ~ 0.6, while $0 — 0.07 
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TABLE IV. Low-temperature values of the surface tension at 
T = 2MeV (in practice, similar to T = 0). We choose a 
sigma meson mass m CT = 500 MeV and pure glue transition 
temperature To = 182 MeV. The parameters not shown in 
this table are found in Table HI 



Parametrization 


Log [19] 


Poly-I [17] 


Poly-II [18] 


E [MeV/ftn*] 


13.0 


13.0 


28.2 



for the Poly-II parametrization and $o — 0.5 for the 
Log parametrization. A comparison to the pure gauge 
SJ7(3)-lattice calculation of [76|, [TtJ shows that one has 

The reason why our results for the surface tension with 
the Poly-II parametrization arc so different can be un- 
derstood from the definition (|3.21|) . The kinetic param- 
eter k in (|3.15j) is fitted from the value of the surface 
tension in pure gauge ST/(3)Tattice simulations using 
Eq. p,14[) . After some simple manipulations, one sees 
that the smaller the barrier between the minima of the 
pure gauge potential, the larger k has to be in order to 
reproduce the known value of S S c/(3) ^ 0.016T 3 [z34z3 - 
Once the Poly-II parametrization leads to a very short 
and narrow barrier, it gives the large value n ~ 160 To. 
This is three orders of magnitude larger than the values 
given by the Poly-I and Log parametrizations. As a re- 
sult, a large value of k leads to a large coefficient h (|3.22p 
for the Poly-II parametrization and, consequently, to the 
behavior of the surface tension seen in the inset of Fig. [5] 
In the following, we only discuss the surface tension de- 
rived from the Poly-I and Log parametrizations. 

Our results on the zero-temperature surface tension 
given in Table IIVI and the temperature dependence of 
the surface tension seen in Fig.[B]are very similar to those 
found in Refs. [3(1 HH, which considered the two-flavor 
QM and NJL models, respectively. This means that the 
addition of the strange sector does not change apprecia- 
bly the dynamics of the phase transition at low temper- 
atures and high chemical potentials. Notice that, given 
a point of phase coexistence in the T — [i plane, Aa x is 
much larger than Aa y , except for conditions very close 
to the CEP. Therefore, the coefficient h in Eq. at 
low temperatures is barely changed with the contribution 
from the strange sector. Actually, at low temperatures, 
the Polyakov loop sector does not give any sensible con- 
tribution to the surface tension either. The temperature 
range in Fig. [5] where the profiles of the surface tension 
with the Log and Poly-I parametrizations differ indicates 
where Polyakov loop sector contributes to the surface ten- 
sion and the difference of the profiles can be considered 
as an estimate of the uncertainty of the Polyakov loop 
contribution. 



A. Implications for proto-neutron stars, the early 
Universe and heavy ion collisions 

The values of surface tension for the Nf = 2 + 1 PQM 
model wc found have interesting implications for several 
physical scenarios. For example, compact stars can be 
considered as laboratories for nuclear matter at low tem- 
peratures and at such high densities that they may con- 
tain quark matter [78||. Possible scenarios for the forma- 
tion of quark matter in compact stars are old accreting 
neutron stars, proto-neutron stars after a supernova ex- 
plosion or during the early postbounce evolution of core 
collapse supernovae [79l |. Physical conditions and time 
scales in these cases imply equilibrium with respect to 
weak interactions and low electron fractions. Estimates 
from [27l | show that a hadron-quark phase transition dur- 
ing the bounce phase of a core-collapse supernova can 
be dynamically suppressed if the surface tension of this 
phase interface is much larg er than, say, 20 McV/fm 2 . 
The estimates from [10, [3lf and this work consistently 
point towards low values of the surface tension, which 
would be compatible with the formation of quark mat- 
ter during the bounce. An observable signal would be 
a second peak in the neutrino signal dominated by the 
emission of antineutrinos and with a significant change 
in the energy of emitted neutrinos [79l |. However, none 
of these calculations really takes into account realistic 
equations of state for supernova matter, which has to 
include not only scalar mesons, but also vector mesons, 
nucleons and, very importantly, leptons. A calculation of 
surface tension in such a complete model would be most 
welcome. 

In the cosmological case, physical boundary conditions 
to describe the QCD phase transition in the early Uni- 
verse include charge neutrality, equilibrium with respect 
to weak interactions and baryon and lepton asymmetries 
consistent with observations. Observations of the cosmic 
microwave background radiation and constraints from 
primordial nucleosynthesis require a tiny baryon asym- 
metry at relatively low temperatures (T < 1 MeV) [8(| . 
In the standard scenario this observational constraint is 
extrapolated up to the scale of the QCD phase transition 
that is then a smooth crossover. In the scenario of lit- 
tle inflation [28|, the Universe enters the QCD era with 
a very high quark density. As a result, the quark-gluon 
plasma (QGP) that fills the Universe cools down until it 
eventually crosses a first-order line of the phase diagram, 
becoming metastablc. While the QGP is mctastable, 
the conditions for a cosmic inflation can be met and 
the extra quark density becomes very diluted. Obser- 
vational signals of this evolution include an enhancement 
of primordial density fluctuations on stellar up to galac- 
tic scales, production of galactic and extragalactic mag- 
netic fields and a modification of the gravitational wave 
spectrum 0, [8l[ . At some point of the expansion, how- 
ever, the phase transition from a QGP to hadronic matter 
must happen, most likely through bubble nucleation. In 
order to be effective, the baryon dilution required by the 
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little inflation scenario needs to be long enough. This 
requires that the QGP remains metastable even for high 
degrees of supercooling, something that requires a very 
low nucleation rate and, therefore, a large value of sur- 
face tension. As we have discussed, however, the values 
of surface tension found with chiral models, including 
the study we have performed in this work, are relatively 
small and possibly do not allow the strong metastabil- 
ity required by the scenario of little inflation. However, 
also a large lepton asymmetry can drive the evolution- 
ary path of the Universe towards larger quark chemical 
potentials (S2| | . 

In heavy ion collisions the phase boundary of hadronic 
and quark matter is, if at all, crossed twice. First, as 
formation of a quark gluon plasma in a hadronic gas and 
then as rehadronization of the fireball. Complications in 
studying nucleation in heavy ion collisions are the short 
time scales and the finite size of the system. Addition- 
ally, the nucleation rate has to be considered in relation 
to the expansion time. These conditions can lead to the 
fact that the system stays in the metastable state close 
to the spinodal instability and that the dominant mech- 
anism for phase conversion is the alternative scenario to 
homogenous nucleation, namely spinodal decomposition 
[23j. Nevertheless, the growth rate of fluctuations by 
spinodal instabilities is closely related to the surface ten- 
sion [83|, HH ■ They can lead to observabl e sig natures for 
the value of the surface tension we found [84|] ■ Addition- 
ally, these fluctuations can be amplified by nucleation 
in the metastable region. The relatively small values of 
the surface tension we found suggest a early nucleation 
of small quark-gluon plasma droplets at relatively mod- 
est energies like at FAIR's SIS 100 [Hf. The details of 
rehadronization leave their fingerprints on those observ- 
ables that are sensitive to the life-time of the fireball. 
Weak supercooling favors the thermal freeze-out to hap- 
pen in the hadronic phase with impact on particle yields 
and spectra (86j and a distinct hydrodynamic expansion 
pattern [23j . 



V. CONCLUSION 

In this work, we have analyzed some thermodynamical 
properties of the PQM model with Nj = 2+1 constituent 
quarks. In particular, we were concerned about the defi- 
nition of the equilibrium state of a system at temperature 
T and quark chemical potential fi. As in the case of the 
PNJL model, the in-medium effective potential of the 
PQM model is not a real function of real variables and, 
therefore, it has no minima. As a way to circumvent this 
problem, we have first rewritten the effective potential in 
terms of real variables only, such that the real and imag- 
inary parts of the potential are separated. The relation 
between the imaginary part of the effective potential and 
the fcrmion sign problem is discussed. The real effective 
potential is found to be the appropriate quantity to be 
minimized and it is a consistent quantity in terms of stan- 



dard arguments of equilibrium statistical physics. This 
allowed us to calculate some properties of the model in 
equilibrium at finite T and \i. In particular, we calculated 
the evolution of the order parameters as functions of the 
temperature (for fj, = 0) and of the quark chemical po- 
tential (for T = 2 MeV, which is, in practice, equivalent 
to T = 0). We also calculated the T — fi phase diagram 
of the model and compared our results to the ones found 
in the literature using the saddle-point method. The re- 
sults are equivalent for fi = (where they had to be), 
but differ for nonzero fi. The phase transition lines of 
both methods coincide but the location of the critical 
end point can be at up to 40% larger chemical potentials 
using the minimizing method depending on the Polyakov 
loop potential parametrization and parameters. 

A careful minimization of the effective potential also al- 
lowed us to study the problem of homogeneous nucleation 
of bubbles in a first-order phase transition. More specif- 
ically, we calculated the surface tension of an interface 
between the two phases predicted by the model, a quan- 
tity that is crucial for the nucleation rate of bubbles in a 
first-order phase transition. We saw that the Nf = 2 + 1 
PQM model yields results very similar to those of the 
two-flavor NJL and QM models, so that the influence 
of both the strange quark and the Polyakov loop at low 
temperatures is small. However, the same cannot be said 
for higher temperatures, where all contributions become 
of the same order of magnitude. Our overestimate gives 
a conservative upper bound of S < 15 MeV/fm 2 for the 
surface tension. The actual value, however, may be even 
lower, not only because of the direct approximations in 
the calculation of the surface tension, but also due to vac- 
uum terms, which we have neglected and should make 
the first-order transition weaker at low T and high /1. 
In any case, this reinforces the trend shown recently in 
Refs. [10, HH of a low surface tension in chiral models for 
QCD at finite baryon density. Such a low value would al- 
low a quick hadron-quark phase conversion. This implies 
interesting implications for several physical scenarios, be 
it heavy ion collisions, proto-neutron stars or the early 
Universe. 

In summary, the analysis carried out in this work sug- 
gests that, in spite of the good agreement of the chiral 
models at finite T and /i = with lattice calculations, one 
should take care when the fi 7^ case is addressed. This 
difficulty is manifest in the non-reality of the equilibrium 
effective potential for /j / 0, even far from any phase 
transition. Notice that, in principle, the same problem 
can affect any chiral model with gauge fields coupled to 
quarks at finite fx, such as the PNJL model. We be- 
lieve that more consistent solutions to these (and other 
possible) problems are still needed as they would bring 
more confidence to further progress in the domain of high 
chemical potentials of the QCD phase diagram with ef- 
fective models. 
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